Generalising the logistic map through the (/-product 



R W S Pessoa, E P Borges 

Escola Politecnica, Universidade Federal da Bahia, Rua Aristides Novis 2, Salvador, Bahia, 
Brazil 

E-mail: robsonpessoa2007@gmail . com, ernesto@ufba.br 

Abstract. We investigate a generalisation of the logistic map as x n +i = 1 — ax n ®q map x n 
(—1 < x n < 1, < a < 2) where ® g stands for a generalisation of the ordinary product, known 
as q-product [Borges, E.P. Physica A 340, 95 (2004)]. The usual product, and consequently 
the usual logistic map, is recovered in the limit q — > 1, The tent map is also a particular case 
for QWp — ¥ oo. The generalisation of this (and others) algebraic operator has been widely used 
within nonextensive statistical mechanics context (see C. Tsallis, Introduction to Nonextensive 
Statistical Mechanics, Springer, NY, 2009). We focus the analysis for q ma p > 1 at the edge of 
chaos, particularly at the first critical point a c , that depends on the value of q m a P - Bifurcation 
diagrams, sensitivity to initial conditions, fractal dimension and rate of entropy growth are 
evaluated at a c (q map ), and connections with nonextensive statistical mechanics are explored. 



1. Introduction 

Low-dimensional non-linear maps represent paradigmatic models in the analysis of dynamic 
systems. The discrete time evolution and the small number of relatively simple equations make 
their treatment easy, without losing the richness of the behaviour, exhibiting order, chaos and 
a well defined transition between them (see, for example, [H[2])- 

Strongly chaotic systems are of special interest for statistical mechanics, once they feature well 
known characteristics: exponential sensitivity to the initial conditions, ergodicity, exponential 
relaxation to the equilibrium state, gaussian distributions [3]. 

In-between ordered systems (with negative Lyapunov exponent) and (strongly) chaotic 
systems (with positive Lyapunov exponent) there are those with zero maximal Lyapunov 
exponent. These systems are characterised by power-law sensitivity to initial conditions, instead 
of the exponential sensitivity, and thus are considered as weak chaotic systems. This change in 
the dynamics may lead to break of ergodicity, non-exponential relaxation to equilibrium and/or 
non-gaussian distributions. These behaviours are usually expected to be found in systems that 
are described by nonextensive statistical mechanics [HE]. Some low dimensional maps, e.g. the 
logistic map, also exhibit weak chaoticity at the edge of chaos, and hence the interest in studying 
them to better understand nonextensivity. 

Power-law like sensitivity to initial conditions and power-law like relaxation to the attractor 
(more precisely a ^-exponential law) have already been found in logistic-like maps [21 [6]. q- 
exponential function {e x q = [1 + (1 — q)x] 1 /^ 1 9 , the subscript + is explained in the following) 
appear within nonextensive statistical mechanics and it generalises the usual exponential 
function (recovered as q — > 1). It is asymptotically a power-law (for q > 1 and x < or 
q < 1 and x > 0). Sensitivity to initial conditions of the logistic map at the edge of chaos is 



identified to a (/-exponential, with a specific value of the parameter q, denoted as Qserf The rate 
of entropy growth S q /t (S q is the nonextensive entropy, defined later by Eq. Q, and t is time) 
is another parameter usually evaluated in maps. It must be finite at the macroscopic limit, and 
there is one special value of q denoted q en t (from entropy) that makes S qent /t finite. At the edge 
of chaos, q en t ^ 1. It is numerically verified that q sen = q en t (see [5] and references therein). 
Relaxation of the logistic map to the attractor at the edge of chaos also follows a (/-exponential 
behaviour, with a different and specific value of the parameter q, denoted q re i. The relation 
between q en t < 1 and q re i > 1 plays a central role in the foundations of nonextensive statistical 
mechanics. For completely chaotic systems, these values collapse to q en t = q r ei = 

1 (See [5] for 

details) . 

Nonextensive statistical mechanics has lead to developments in many related areas, including 
generalised algebras [3 [8]. These works have introduced generalised algebraic operators, and 
here we are particularly interested in the g-produclQ defined as 



i 



x ® q y = sign(x)sign(y) 9 + \y\ l q — 1 1 q (1) 

where the symbol [A]+ means that [A]+ = A if A > and [A]+ = if A < (known as 
cut-off condition, [A]+ = max{0,A}, for short). The limit q — > 1 recovers the usual product 
(x <8>i y = xy). Our work consists in generalising the logistic map as 

x n+l = 1 - a (x n ® qmap x n ) ^ 

= 1 - o[2|x„| 1 -*»«f - 1] + 9map 

(-1 < x n < 1, < a < 2). The cut-off condition implies that [2\x n \ l ~ qmav - 1]+ = 0, and 
thus x n +i = 1, if \x n \ < \/2 l K l ~ qmap \ This (/-logistic map recovers the usual logistic map 
for q map = 1, and also the tent map for q m ap — > +oo (the tent map properly shifted as 
x n+ \ = 1 — a|x n |). At the limit q map — > — oo, it becomes x n +i = 1 for — 1 < x n < 1 and 



x n +i = 1 — a for x n = ±1. Some details regarding these limits are sketched in Appendix A 
Figure [T] shows one iteration of the map. 

The z-logistic map (x n+ \ = 1 — a\x n \ z , z > 1, < a < 2) [TU [TJl E31 E], is another 
generalisation of the logistic map for a general power z > 1 and holds some similarity to the 
present g-logistic map. 

Bifurcation diag rams for different values of Qmap shown in Fig. [2j As q map goes from 1 (the 
usual logistic map) to 2, the value of the parameter a for the first bifurcation goes from a = 0.75 
to a = 1 (it approaches 1 from the left). Also the value of the parameter for the accumulation 
of bifurcations goes from a c = 1.401155189092. . . to a c = 1 (it approaches 1 from the right). 
It means that the period doubling cascade becomes narrower until it eventually disappears at 
Qmap = 2. Similar behaviour also happens with the other windows of order inside chaos (where 
there is tangent bifurcation): they get narrower as q ma p increases, and eventually disappear for 
Qmap = 2. Complete chaos is preserved at a = 2, Vg map . The Schwarzian derivative for the 



g-logistic map is given by (see Appendix B ) 



(57)0*0 - [Qmap 2 q m ap J 4 _ 4 | x | gmap -l + \ x \2(q map -l) 

This expression is negative in the interval < q ma p < 2, and positive for q map > 2 and for 
qmap < (SD(f(x)) = for 

qmap — and for q m ap — 2). This means that the route to chaos 
for < q m ap < 2 is by period doubling bifurcation. 

1 The g-product has been also used in the generalisation of Gauss's law of errors [9], in the formulation of the 
q- Fourier transform, and in the generalisation of the central limit theorem [10] , 
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Figure 1. x n +i as a function of x n for the 
g- logistic map (Eq. ^ with a = 2). The 
usual parabolic behaviour is recovered at 
Qmap = 1- Tent map is found at q — > oo 
(with the peak parameter set to unit). As 
Qmap departs from 1 towards — oo (q m ap < 
1) the cut-off condition in Eq. Q yields an 
increasing region in which x n +i = 1. For 
Qmap — oo, x n> o alternates between 1 
and (1 - a) Vx £ [-1,1]. For q map > 1 
the map is discontinuous at x = and 
for qmap < 1 the discontinuity is at |x| = 
l/2i/(i-9). 
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Figure 2. Bifurcation diagrams for differ- 
ent values of q ma p (indicated). Windows of 
order inside chaos vanishes as q ma p — > 2. 
In this paper we explore q map > 1 but 
Qmap = 0.5 is shown as an instance just to 
give an idea of the scenario for q map < 1: 
the regions of chaos become narrower and 
the regions of order become dominant. 



2. Sensitivity to initial conditions 

The maximal Lyapunov exponent may be evaluated at the edge of chaos according to (see, for 
instance, [3]) 

-. N-l 

X mca = lim -£ln !/'(*,)! (4) 

iV— >oo iV n 

where f'(x) is the derivative of the q- logistic map, 

Qmap 

f(x) = -\x\'2a\x\- qma v 2|x| 1 ~ 9map - 1 ^ qmap . (5) 

\x\' is the derivative of the absolute value function: \x\' = 1 for x > and \x\' = —1 for 
x < 0. For values of the control parameter a different from critical ones, the sensitivity to initial 
conditions are characterised by exponential divergence at regions of chaos and exponential decay 
at regions of order, i.e., positive or negative Lyapunov exponent Ai (subscript 1 will be clear in 
the following) in 

m= lim ^ t \ = e^ t . (6) 
sv 1 Ax(oHo Ax(0) v ; 

At the edge of chaos it was proposed that the divergence follows an asymptotic power-law (in 
fact a ^-exponential law) [2] characterising a slow dynamics, 



(7) 
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Figure 3. Sensitivity to initial conditions as defined by variable L (Eq. ([8])) for different values 
of q m ap (indicated). 40 initial conditions are used. For q ma p = 2.5 (figures e and f) two values 
of the control parameter a are presented. It can be seen that it is possible to have coexistence 
of attractors for a certain ensemble of initial conditions (for a = 1.001), one with L > and 
the other with L < 0. The same value of q ma p = 2.5 and with a = 1.1, the initial conditions 
always lead to L > 0. Fig. f shows linear dependence of L with N . Slopes (in absolute values) 
for a = 1.001 for L > and L < differ. 



where Ax(0) represents the distance between two neighbouring initial conditions and sen stands 
for sensitivity to initial conditions (q se n < !)• Eq. Q is recovered at q sen = 1 (see Appendix 



[A| and this is the reason for the subscript 1 in Ai, Eq. ([6]). As a graphical representation of the 
sensitivity, it can be defined the variable L as [2]: 



N-l 



L=J2^\f(xi)\- (8) 



i=0 



Figure [3 shows five instances of L vs. IniV. For q map = 2 (Fig. [3ji) the dependence of L on N 



is very s 



. ow and cannot be seen up to ./V = 2 15 (the upper limit of the figures). For q map > 2 it 
is possible to have coexistence of attractors according to the control parameter a: depending on 
the initial conditions, L may be increasingly positive or decreasingly negative. Fig[3f shows that 
there is a linear dependence of L on but with different slopes for the cases L > and L < 0. 
For values of a > 1.07499 . . ., L < is never exhibited. Coexistence of attractors was also found 
in |15] for a different deformation of the logistic map (the authors also call their deformation as 
q- logistic map). 

Lyapunov exponents are displayed in Fig. [4] showing transitions from order to chaos. The 
figures show that these transitions become sparse as q ma p departures from unit, and for q m ap = 2 
there is only one transition (robust chaos) [16] . Fig. [4]f shows coexistence of attractors for 
qmap = 2.5. 

Figure [5] illustrates the first point of accumulation of period doubling bifurcation 
function of q ma p- For 1 < q map < 2, the behaviour is ordinary in the sense that for a < a c , 
L < 0, and for a slightly greater than a c , L > 0. For q ma p > 2 a different behaviour appears: 
order is found for a < 1 (region below the solid line), while chaos is found above the dashed 
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Figure 4. Lyapunov exponent as a function of the control parameter a for different values of 
Qmap (indicated). Up to the time limit used in Eq. Q, the almost vertical line in Fig. c at 
a = 1.7477 does not crosses zero, but it does present negative values if the calculation is done 
with a greater precision (it was used t = 2 18 and Aa = 2 x 10 -4 ). Fig. f exhibits coexistence 
of attractors for two different initial conditions. It's inset is an amplification of the region of 
coexistence of attractors for two initial conditions (indicated as full circles and open squares). 
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Figure 5. Dependence of the first point of accumulation of period doubling bifurcation a c on 
Qmap- For q m ap > 2 there is a region that exhibits coexistence of attractors depending on the 
initial conditions. Values of L for this figure were calculated without transient time and final 
time t en d = 2 18 , and increments for the control parameter Aa = 10~ 7 . Critical point is taken as 
that for which < A moT < 5 x 1CT 5 . 



line. The region in-between presents coexistence of ordered and chaotic behaviours, depending 
on the initial conditions (see figures [3^, [3^ and|4]t). 

Table [T] shows the range of the critical points a c (first point of accumulation of bifurcations) 



Table 1. Critical points 
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for different values of q ma p- The value of a c is between o_ and a+. The fourth column shows 
adopted value for a c . For the evaluation of a c listed in Table [TJ the Lyapunov exponents were 
calculated with a transient time of ttrans = 2 23 and a final time of t en d = 2 23 (the values of the 
Table are more accurate than those of Fig. [5]). Initial condition was fixed in xq = 0.65 for all 
cases. 

In Fig. [6] we show the cycle 3 window for q ma p = 1.25. This window of order inside chaos 
(as well as all the others) becomes narrow: in order to identify the Lyapunov exponent in this 
instance it was necessary to give increments of Aa = 1.0 x 10 -8 and t en d = 2 23 , with a transient 
ttrans = 2 23 . Identification of windows of order in this q- logistic map is computationally time 
consuming as q map departs from unit. 

3. Entropy production 

Parameter q en t (from entropy) in Tsallis entropy [4J 

q-1 

(W is the number of microstates and pi is the probability of microstat^] i; we use k = 1 
without loss of generality) is estimated according to the method developed in [T7] . It consists 
in calculating the rate of increase of entropy, which must be finite for a large (virtually infinite) 
system^] The phase space (x n 6 [—1, 1]) is divided into W cells, with ./V points (initial conditions) 
inside one of them. As time evolves, Tsallis entropy S q is calculated for many values of q. 
Calculation is done again with N points inside another initial cell, and this process is repeated 
for each cell of a certain ensemble called "best initial condition cells" (the definition of "best 
initial condition cells" is: the integrated number of occupied cells must be greater than a certain 
(arbitrary) threshold, the most visited cells). Then it is taken the average of entropy for each 
time step, and this average is finally plotted against time, for various values of q (see Fig. [7] 

2 In statistical mechanics, for a given macrostate (that is, macroscopically measured state) there are a number 
W of compatible configurations of the constituent elements, or microscopic states (microstates, for short). A 
macroscopic variable is, in fact, a measure of the average of all compatible microscopic states. 

3 By large system we mean a macroscopic system. The density variable of a quantity F, F/N, must remain finite 
as the number of microscopic constituents TV reaches the thermodynamical limit limjv^oo- 
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Figure 6. Maximal Lyapunov exponent for q ma p = 1-25 at the vicinity of the largest window, 
of cycle 3, with increments of Aa = 1.0 x 10~ 8 . It was used a transient time of ttrans = 2 23 
and the time used to estimate the Lyapunov exponent t en d = 2 23 (with lower time intervals 
the transition from negative to positive Lyapunov exponent cannot be seen). Inset corresponds 
to the dashed rectangle in the main panel and shows points of zero Lyapunov exponents that 
corresponds to period doubling bifurcation. 




for an instance). There is only one special value of q for which the increase of S q is linear in 
the macroscopic limit, that is, the production of entropy n g = lim^oo limjy->oo um iV->cx> &q/t 
remains finite. Xhis special value is identified with q en t (see P2] for details; in that paper, q ent is 
called q*). For values of the control parameter a that corresponds to chaotic behaviour q en t = 1, 
so ordinary Boltzmann-Gibbs-Shannon entropy is the proper one to be used. This is not the 
case at the edge of chaos (a = a c ), which we are interested in this work, and the value of q en t 
which yields linear growth of entropy (and consequently finite rate at the macroscopic limit) is 
smaller than one (for the ordinary logistic map, q ma p = 1, at the edge of chaos, q en t = 0.2445...). 
We applied this procedure to the (/-logistic map. The result of this procedure for various values 
of q m ap is presented in Fig. [8j 

The integrated number of occupied cells, used to identify the best initial condition cells, is 
shown in Fig. [9j We see a reduction of the number of visited cells as q ma p goes from one to two 



and a change in the displayed pattern. Fig. 10 shows number of cells with visits greater than 
5000 as a function of q ma p- 

4. Relaxation to the critical attractor 

Relaxation to the critical attractor was presented at [6] and basically consists in the division 
of the phase space [-1,1] into W cells and take an ensemble of N initial conditions uniformly 
distributed in the entire phase space (thus maximal entropy). Time evolution leads to the 
decreasing of the number of occupied cells W OC c(t) according to a power-law with log-periodic 
oscillations (W occ (0) = W). It is supposed |6] that the power-law is the asymptotic limit of a 
^-exponential, and the parameter is denoted as q re i, for relaxation (q re i > 1): 



W occ (t) = (1 + (1 - q r el)K qr j) ^rel , (10) 



K q is the inverse of a characteristic time. Figure 11 presents the results for the fraction 



of occupied cells W occ (t) /W occ (0) for different values of q ma p at their corresponding critical 




Figure 7. Rate of increase of entropy for q ma p = 1.1. We found q en t = 0.33. This value is found 
by fitting a parabola (S q (t)) = a + bt + ct 2 for the time interval 15 < t < 38 (the same used in 
[T7]). Initial times are always left out of calculation because the concavity is always positive at 
that region. For q < q ent the concavity of the curve is positive (squares, c > 0) and thus the 
rate of growth of entropy diverges at the macroscopic limit, which is unphysical. For q > q en t, 
the concavity is negative (triangles, c < 0), thus the rate of growth of entropy is zero at the 
macroscopic limit, which is also unphysical. c = is shown as circles. 
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Figure 8. q en t as a function of q ma p- Inset shows abscissa in a different scale (note that abscissa 
is in log scale). The value at q ma p = 1-02 may have numerical imprecision. Lines are only guide 
to the eyes. 



points. Log-periodic oscillations present increasing periods for q m ap — > 2. The slope in the 
log-log plot (Fig. |IT^ ) at the region of the log-periodic oscillations is used to estimate qrel &s 
slope = 1/(1 — q re i)- Fig. 11 3 presents q re i as a function of q ma p- 

The procedure for estimating q re [, that is, the shrinking of the number of occupied cells, 
is also used to estimate the fractal dimension df at the edge of chaos (it is a kind of box 
counting method). JVo CC (oo) increases with W occ (0) according to W occ (oo) oc [W O cc(0)] a! ' f • Fractal 



dimension decreases with q ma p as shown in Fig. 12 
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Figure 9. Integrated number of vis- 
its per cell for five cases {q ma p = 
1, 1.30, 1.50, 1.70, 1.90). Phase space is di- 
vided into W = 10 5 cells, each one contains 
10 6 points uniformly distributed. Time 
evolves up to 50 iterations. Figures show 
x 6 [0.5,0.6] just for better visualisation. 



Figure 10. Number of cells with more 
than 5000 visits in a run of 50 iterations 
as a function of q map - N pc = 10W is 
the number of initial conditions per cell. 
Dashed horizontal line in Fig. [9] indicates 
the threshold of 5000 visits per cell. 
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Figure 11. Relaxation to the critical attractor at the edge of chaos. Fig. a is a log- log plot of 
W occ (t) /W occ (0) x t for three different values of q map . It can be seen the results of regression of 
power-laws, whose slopes are identified as 1/(1 — q re i)- Fig- b shows q re i as a function of q ma p- 
It was used W = 2 14 and number of initial conditions per cell N pc = 2 13 . Lines are only guide 
to the eyes. Fluctuations in the curve of Fig. b are due to inaccuracy — the results are very 
sensitive to small changes in the characteristic parameters. 



5. Final remarks 

The generalisation of the logistic map by means of the (/-product introduces some interesting 
features in its dynamical behaviour. Our analysis is focused on q ma p > 1. In this region, the 
windows of order inside chaos become narrower as q ma p increases until all of them disappear 
and the map becomes the tent map. For q ma p < 1 (not analysed in this paper) the opposite 
behaviour occurs: the regions of chaos become narrower and order dominates the scenario. A 
remarkable feature of the (/-logistic map is to continuously pass from a map with a variety of 
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Figure 12. Fractal dimension at the edge of chaos as a function of q ma p. For the usual logistic 
map {q m ap = 1), df{a c ) = 0.53665. Lines are only guide to the eyes. Fluctuations in the curve 
are due to inaccuracy — the method is very sensitive to small changes in the characteristic 
parameters. 



behaviours, such as period doubling bifurcation, multifractality and power-law like sensibility to 
initial conditions at the edge of chaos (the logistic map), to a robust map (the tent map). We 
have calculated the sensitivity to initial conditions, the rate of entropy growth, the relaxation 
to the critical attractor and the fractal dimension at the edge of chaos. These methods permit 
to estimate the parameters q en t and q re i for different values of q m ap- The entropy parameter 
and the relaxation parameter are two indices that appear within the context of nonextensive 
statistical mechanics and the understanding of their dependence on the control parameters of 
the system may lead to their a priori determination. Some other evaluations for this (/-logistic 
map remains to be done, e.g. the dependence of q re i on the coarse graining W and its relation 
to q sen (as in |12] for the z-logistic map, and as in |18l I19j for the Henon map), the probability 
distributions of sums of iterates as in |20|I21|. multifractality as it was done in [22J, and tangent 
bifurcations. 
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Appendix A. Special limits for the g-product 

The limit q — >• v where v is one of the following {— oo, l,oo} for the g-product x ® q x as it 
appears in the g-logistic map, Eq. ^ and also in Eq. ([!]) with y = x, leads to {0 ,l°°,oo } 
respectively. The indeterminates can easily be solved by means of a simple trick that is to 
rewrite the g-product as 

x® q x = [2\x\ l - q - l)f~ q =h(q) 9 ^ (A.l) 

with h(q) = [2|x| 1-9 — 1] + (note that h(q) > due to the cut-off condition in Eq. 0) and 
9(q) = j^j- Of course h is a function of q and x but for now we are not interested in the 



dependency on x, so we omit it. We also omit the subscript of q for the sake of brevity. Then 



i 



lim^lxl 1 -" - 1]|" 9 = e iim 9 ^9(g)inM<?). (A.2) 

Straightforward application of L'Hospital rule leads to (we remind the reader that \x\ < 1) 

1, if \x\ = 1 1 r 

if \x\ < 1 f fOT V = 
lim x® a x = I o r i (A. 3) 

q^v y ] x z , for V = 1, 

for w = co. 

L'Hospital rule must be applied twice in the case v = oo. A similar procedure applied to Eq. Q, 
now with = 1 + (1 - q S en)\ se J and g(q) = 1 _^ en , leads to Eq. 

Appendix B. Schwarzian derivative for the (/-logistic map 

The Schwarzian derivative is defined by 

(Sf)(x) = l^-U^j . (B.l) 

The function f(x) = 1 — a (x ®q map x) represents the g-logistic map and it may be written as 
f(x) = 1 — ah(x) 9 ^ with h(x) and g(q) given by Eq. (A.l) — now we are interested in the 
dependency of h on x. The three first derivatives of f(x) are given by 

fix) = -a^h&ti 



-aj^Sj&^h'h" 
-cijt^h^-ih 



with h' = |z|'2(l - q)\x\-*, h" = -2q(l - ^)|£c|~ 1 - £? and ti" = \x\'2q(l - q)(l + q)\x\- 2 -i. 
Substitution in Eq. (B.l) leads to 



| 2 ^[2|x| 1 -'? — |;r| 1 +9[2|x| 1 - ( ? - 1]+ |x| 2 

2q c ^ 2 



2 VM^Ix! 1 -? - 1] + |x| 



(B.2) 



This expression may be rearranged as in Eq. ^ that is more convenient to analyse its signal. 
The Schwarzian derivative for the g-logistic map presents a divergence at x = for 1 < q < 2 
and a divergence at x = and at x = ±l/2 1 /i 1 ~^ for q < 1. The case g ma p = 1 corresponds 
to the Schwarzian derivative for the usual logistic map, (Sf){x) = — 3/(2x 2 ), as can be easily 
verified. 
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